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The distributions of the quark number and chiral condensate over the gauge fields are computed for 
QCD in Euclidean space at nonzero quark chemical potential. As both operators are non-hermitian 
the distributions are in the complex plane. Moreover, because of the sign problem, the distributions 
are not real and positive. The computations are carried out within leading order chiral perturbation 
theory and give a direct insight into the delicate cancellations that take place in contributions to 
the total baryon number and the chiral condensate. 
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I. INTRODUCTION 



The phase diagram of strongly interacting matter is determined by the behavior of the chiral condensate, (a), and 
the quark number, (n). These quantities signal the breakdown of chiral symmetry and the formation of baryonic 
matter. Not only the expectation values of a and n are of interest, but also their fluctuations are central to our 
understanding of strongly interacting matter since they may be visible in fluctuations of the observables measured in 
heavy ion collisions [H, U d, Q ■ 

In Euclidean space both a and n are complex when evaluated at nonzero chemical potential for a typical gauge 
field background. The reason why n and a take complex values is that the fermion determinant at nonzero chemical 
potential is complex: if we write det(D + fjryo + m) = rcxp(i9) we have 

d . d d 

n = — log dekiD + ^70 + m) = — log r + i—0 (1) 
d/j, d/i d/j, 

d . d d 

o = - — logdet(£> + /J7o + m) = —— logr + 1—— 0. 
dm dm dm 

The fluctuations, therefore, take place in the complex a- and n-plane. Since the fermion determinant is also complex 
valued for non zero chemical potential, fluctuations in all directions of the complex a- and n-plane contribute to the 
expectation values (a) and (n), which of course are real. A description of fluctuations of a and n is contained in 
their distributions, P„ and P n . Here we derive these distributions analytically to leading nontrivial order in chiral 
perturbation theory. From ([T]) we see that the imaginary part of the quark number operator is directly related to 
the phase of the fermion determinant. We therefore expect that the cancellations due to fluctuations of the phase 
are tightly linked to the distribution of the quark number in the imaginary direction. The results presented below 
confirm this expectation. 

Chiral perturbation theory Q is the low energy limit of QCD which describes the strongly interacting theory 
in terms of weakly interacting Goldstone modes corresponding to the spontaneous breakdown of chiral symmetry. 
Since the pions have zero quark charge, the expectation value of the quark number in chiral perturbation theory is 
automatically zero. The square of the quark number 

is of course also zero since in chiral perturbation theory the partition function, Z , is independent of fx. The distribution 
of the values which the quark number operator takes as the gauge fields fluctuate is, however, nontrivial even when 
evaluated within chiral perturbation theory. To see that this is necessarily the case, let us consider the second moment 
of n i.e. (n 2 ) = j dxdy(x + iy) P n (x + iy). If the second moment is non zero the distribution P n is nontrivial, i.e. not 
a two dimensional 5-function at the origin of the complex n-plane. 

Let us emphasize that (n 2 ) is not the average of the square of the quark number, see (J5J), rather the second moment 
must be thought of as 

Z dfj, u d/j, d 

evaluated at degenerate chemical potentials and quark masses for the up and down quark. Even though there is an 
equal number of quarks and anti-quarks associated with each pion, the vev of (n 2 ), is nonzero if there are correlations 
between the quarks and the anti-quarks of different flavor. Confinement of quarks and anti-quarks into pions strongly 
suggests that such a correlation exists in chiral perturbation theory, and in fact (n 2 ) takes a non zero value even at 
[i = 0, sec for example Q. At zero chemical potential (n 2 ) coincides with the off diagonal quark number susceptibility 
which has been computed analytically by high temperature perturbation theory [7J as well as by numerical lattice 
simulations [![. 

Since this is a likely source of confusion let us again stress that P n gives the distribution of n, defined in ([T]), over 
the ensemble of gauge fields: If one makes a frequency plot of the values obtained for n in a lattice QCD ensemble 
of gauge fields the shape which emerges is described by P n . The first moment measured on this distribution gives 
the average quark number. The second moment measured on this distribution, however, only give the first of the two 
terms which combine to give the square of the quark number, cf. @. If we phrase this in terms of the eigenvalues Zk 
of 70 (D + m) we have 



(4) 
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So (n^) is not the average of a square and consequently not the second moment of a distribution over the gauge fields. 
The reason is that r? q contains correlations due to the Pauli principle. 

On a more technical level, the reason for the /x dependence of P n (x + iy) can bee seen from the generating function 
of the distributions. These involve quarks with different chemical potentials, i.e. the isospin chemical potential, 
Mi — A*2j is nonzero and couples non-trivially in the generating functions. This point will be explained explicitly 
below. The success of chiral perturbation theory in predicting the behavior of QCD at nonzero chemic al p otential 
has been demonstrated clearly by recent comparisons to lattice QCD results both in the p- regime [H, [HI EI HH and 
the e-regime [H M M, M El • 

The phase of the fermion determinant may lead to exponentially large cancellations in the computation of expecta- 
tion values. If we would try to measure this expectation value numerically with lattice QCD simulations, we would be 
confronted with numerical errors which are exponentially hard to handle. This is the QCD sign problem and severely 
limits first principle studies of the QCD phase diagram (see [l8| for recent reviews of the QCD sign problem). 

The results for the distributions that are derived below give a direct insight in the cancellations caused by the sign 
problem. For example, the vanishing value of the baryon number (within chiral perturbation theory) is obtained only 
after a delicate cancellation between the contribution from Re[n] and from Im[n]. The range of the fluctuations in 
the complex quark number plane which must be taken into account, in order to obtain the total baryon number, 
grows like the four volume. In contrast, the width of the distribution for the ensemble where the phase of the fermion 
determinant is ignored (the phase quenched ensemble) only scales like the square root of the volume. The numerical 
challenge faced by lattice QCD is to go from the phase quenched (simulation) ensemble to the full (target) ensemble. 
One therefore needs to sample the far tails of the distributions. The analytic insight we have obtained here can in 
this way help to understand the limitations of present lattice gauge simulations at non zero chemical potential and 
give hints for future developments. Besides the reweighting method fl^ . l2(il. |2~J] discussed here, the results obtained 
are also relevant for the Complex Langevin method [l^, HH, HH [25|, [26[ as well as for the density of states method 

Throughout the paper we will work with two flavors. We refer to the this theory as the Nf = 1 + 1 theory. The 
phase quenched theory is referred to as 1 + 1* since the absolute square of the fermion determinant corresponds to a 
quark and a conjugate quark. 

This paper is organized as follows. First we consider the region where 2p is less than the pion mass and the mean 
fields do not depend on the chemical potential. We derive the distribution of the quark number operator as well as 
that of the chiral condensate. Results are given both for QCD with dynamical quarks and for the partially quenched 
case. Implications for lattice QCD are discussed in section I VII Before concluding we discuss the distribution of the 
quark number operator and chiral condensate for larger values of the chemical potential. The distribution of the 
quark number for imaginary chemical potential is evaluated in the Appendix. 

II. GENERAL DEFINITIONS AND KNOWN RESULTS FROM CPT 

Before starting the actual computation of the distributions, in this section we introduce notations and explain why 
chiral perturbation theory can give information about QCD at nonzero quark chemical potential. 

The quark number operator is the logarithmic derivative of the fermion determinant with respect to the quark 
chemical potential fj, = fj,i + p-2 

n M = Tr n i 7 ° i ' ( 5 ) 

D + ^70 + m 

while for a we differentiate with respect to the quark mass (we consider degenerate flavors) 

a(ji) = Tr 1 (6) 
D + /i7o + m 

At low temperatures, the QCD partition function and its low energy limit in the form of a chiral Lagrangian are 
independent of \i until the chemical potential is sufficient to balance the energy required to create a baryon. Since the 
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expectation value of the quark number and the chiral condensate are derivatives of the partition function with respect 
to fx, they are independent of fx when evaluated in chiral perturbation theory. On the contrary, expectation values 
that cannot be written as derivatives of the partition function with respect to fx may still depend on the chemical 
potential. This happens when the generating function for the operator under consideration includes quarks with 
different values of the chemical potential. 

To illustrate this let us compute the expectation value of n 2 , the cross correlation introduced above. In order to 
obtain n 2 we start from the generating function 



Zi + i(n, = (det(L> + /i7o + m)det(D + fx^o + m)) 

and differentiate with respect to the two different chemical potentials 

/ 2\ 1 v d d 

( n )i+i = ~ / r I™ — — 

Z i+1 (fx, fi) mi-^m dfx dfii 



(7) 



(8) 



Since the chemical potentials in the generating function are different, there is a nonzero component of the isospin 
chemical potential, fx — fX\. It is this component that induces a nontrivial chemical potential dependence in chiral 
perturbation theory and hence a nonzero result after differentiation with respect to \i and fi\. Since the free energy 
is an even function of the isospin chemical potential (for degenerate quark masses) it takes two derivatives to obtain 
a nonzero value for fi\ — > fx. 

To one-loop order in chiral perturbation theory we have 



Z 1+1 (ix,fx) 

where the one-loop free energy is (the divergent part of Go is independent of \x and does not contribute) 

G {fx,-fx) = ^^2^ ^2^ cosh ( t } ' 



(9) 



(10) 



Thus 



(n 2 )i + i 



hm —-^G {fx,fx{). 



This quantity is usually referred to as the off diagonal quark number susceptibility and is denoted by Xud- 
Below we will use the notation AGq(/xi, ^2) = Go(/ii,/i2) — Go(0,0) and 



(11) 



'X ad 



\ ltd 



- — AG (^i, -fx) 
d/xi 



AG (fXi,Ll 2 

AG (-^i,/i2 



(12) 



dLlldfX2 

d 2 

dLl\d[X2 



for the fx derivatives of the free energy and 

<5S 7 



X B S 



Xs 



—AGo(fx,fx,m,rfi) 
dm 

^AG (fx,-fx,m,m) 
dm 



AGa(fx,fx,m,m) 
AGq(li, —fx, 77i, m) 



(13) 



dmdm 

d 2 
dmdm 



for the one-loop contributions to the chiral condensate and chiral susceptibility. Also for the chiral susceptibility it is 
understood that we will only consider the off diagonal component. Note that all of these quantities are extensive. For 
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quantities labeled by B we have that \x q = [i and /i/ = 0, and for the quantities labeled with I /j q = and fii = /i. 
The quark mass dependence of Go is through via the Gell-Mann-Oakes-Renner relation, 



E(m + to) 
Y 2 ' 



(14) 



Let us also stress that all quantities with superscript B are independent of the chemical potential (in 1-loop chiral 
perturbation theory). 

In table U we give results for one-loop chiral perturbation theory valid for fi < m^/2. The label PQ refers to the 
result obtained in a partially quenched ensemble. By definition this ensemble is generated with the absolute value 
of the fermion determinant but we keep the [i dependence of n and a as in the ordinary theory We will use each of 
these results to check the distributions of the quark number operator and the chiral condensate. 





£ = 1 + 1 


£ = PQ 


(n)e 







(n 2 )s 


/Cud 


v'i + Xud 


\{°)e\ 


E + 25T, B 


E + <5E 7 + 5T, B 


i- 2 )e 


X f + (E + 2<5E S ) 2 


xf + (E + <5E S + ST, 1 ) 2 



TABLE I: The first and second moment of n and a to one-loop order in chiral perturbation theory for fi < ?n^/2. The absolute 
value of the mean field result for the chiral condensate is denoted by E. In the partially quenched (PQ) ensemble the operators 
of the full theory are evaluated for a gauge field background where the phase of the fermion determinant is ignored. 



III. THE DISTRIBUTION OF THE QUARK NUMBER OPERATOR 

In this section we compute the distribution of the quark number in the 1+1 theory. Since 

n{nf = ( Tr - 70 ) = -Tr - 7 ° = -n(-fx) (15) 

\ D + /x7q + to/ D — /i7o + to 

the quark number operator is in general complex (it is purely imaginary at fi = 0). The fluctuations of the quark 
number thus occur in the complex plane. We first derive the distribution of the real part, (n(n) — n(—/j,))/2, and of 
the imaginary part, (n([i) + n(—[i))/2, of the quark number. Then, finally, we compute the full distribution in the 
complex quark number plane. As we shall see the distribution in the complex plane factorizes into the distribution of 
the real part and the distribution of the imaginary part. 

A. The distribution of the real part of the quark number 

Here we derive the distribution of the real part of the quark number defined by 

P££,(s) = (s(x- \{n{v) «(-^)))) i+i , (16) 
within one- loop chiral perturbation theory for (i < m^/2. 
First we represent the (5-function as an integral 

1+1 



P Re[n](z) =^ dk e (e+'jW'-"!-''^ ] . (17) 



Then we expand the exponential of the trace 

( = £V^fr „)_ n( _ M ))A . (is) 

1+1 'f—f 7! \ ' 1+1 



,i|(„( /1 )-„(- (U ))\ ^ (ik/2y 
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This shows that a probability distribution is determined by its moments which we will compute next. 

In order to compute the trace to the jih power we need to introduce 2j replica quarks (see j3l| for an introduction 
to the replica trick in chiral perturbation theory) 

(n^)-n(-^y) i+i (19) 
= lim duj ■ • • eL . ( TT dct(D + ^ 7o + m) 7li det(L> - ^70 + m) ni \ 

rai— >0 m ■ ■ ■ Tlj ' l 11 / 

J u=l / 1+1 

In one-loop chiral perturbation theory the replicated generating function for the real part is given by (this is where 
the assumption fx < m^/2 enters: for /i > m^/2 the generating function is in a Bose condensed phase, see section 

ED 

/jj det(D + Mi7o + m) ni (D - Ati7o + m) n * \ 

\i=l / 1 + 1 

f i 1 \ 

= exp 2J 2n;n ln (G (/i;,/i m ) + G (-fii,fim)) + ^ 2n;(G (/i;, m) + G (-/i;,/i)) . (20) 

U<TO=1 1 = 1 J 

We now take ■ ■ ■ of the generating function, evaluate it at fii = \i and take the replica limits tii — » 0. Note 
that the term linear in the n^s includes Gq(— fJ>i,(J>), and the derivative with respect to \i\ does not vanish at \ii = /i. 
Therefore terms with even as well as odd values of j contribute when we evaluate the derivative d Ml • • • d^ . at Lii = /i. 
All terms with the same number of pairs from the first sum in (|20|) give the same contribution to the jth moment. 
The combinatorial factor for choosing b pairs out of j is 

(21) 

\2b) 6!2" ' (2L) 
For the moments we thus find (recall the notation l|12p): 

IntO/2) / . \ 

((n(/*)-n(-M)) J ') i+i = £ f 2 M(26-l)!![2( x ^ + xL)] fc (2^r 2b 



y [ 3 ) (2^y- 2b - / du (2 w ) 2h e -" 2 /(x^+xL) 



IntO/2) 

1 Z" 00 

= - / cfo (2jBye-^- , ' I ) a /( x - +x «'«), (22) 

where Int(j/2) is the integer part of j /2. These are the moments of a Gaussian distribution centered at vj. The 
distribution of the real part of the quark number is thus given by 

pi+i (x) = 1 P -(—i?/(x B ud +xL) (23) 

Ro N v / , 5 T~ v ; 

We see that f j (a;) is properly normalized. The expectation value of the real part of the quark number equals the 
isospin number in the phase quenched theory evaluated at /i 

/oo 
dxxP^ ] (x)=v I , (24) 
-OO 

and the average of the square of the real part of the quark number is given by 

((RcM) 2 )i+i = dx x 2 P^ n] (x) = uj + \{ X f ld + xD- (25) 

Note that fluctuations of the quark number in the real direction vanish at (i = 0. The reason is that the width of 
the distribution of the real part of n goes to for \x — > 0. 
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B. The distribution of the imaginary part of the quark number 

In this subsection we derive the distribution of the imaginary part of the quark number defined by 

p £w(y) = (f{v+ i\w») + i , (26) 

within 1-loop chiral perturbation theory for [i < /2. 

As in the previous subsection the distribution is determined by its moments which can be shown by representing 
the 5-function as an integral 

> X + X (,,\ - — f°° Ah *- l V k 



P im[n](y) = ^ J dk e-^ k {e*WW-»»)^ (27) 
and expanding the exponential 

( e K«M+»(-,))) i+i = ±U{^(n^+n(-^y\ . (28) 



3=0 J - \ 1 I 1 + 1 

The jth power of Im[n] can be computed by introducing j fermionic replica quarks and j bosonic replica quarks 



it \ , / \vA v 1 a a / ~fj det{D + + my 
{nin + n(-n) ) = hm ■ ■ ■ d H ( \\ — — — — 

/ 1+1 nj->o m • ■ • rij \ dct(L> — /i;7o + my 



1+1 



(29) 



Note that the ratio of the two determinants makes up the phase factor, det(Z?(^))/ det(Z?(— fj,)) = cxp(2i8(fi)). From 
Eq. (fT]) we see that the phase indeed generates the imaginary part of n after differentiation with respect to the chemical 
potential. 

The replicated generating function within one-loop chiral perturbation theory is given by 

/ A det(D + /i t 7o + m) ni \ 
\fJdet(£-Mi7o + m) n */ i+i 

f J i \ 

= exp ^ 2^n ro (Go(/z;,/x m ) - G (-(n,n m )) + y~]2ni(G (fii, fx) - G (-fn, (j,)) 1 , (30) 

\Kto=1 1=1 J 

where we used that /i < m^/2 so that pion condensates are absent in the generating function. 

We now take ■ ■ ■ d^. of the generating function, evaluate it at = \i and take the replica limit m — > 0. The 
only difference with the previous section is the minus sign in front of Gq(~ fii,/i m ) and Go(— /^z, /x). We thus find the 
moments 



Int072) 



((n( M ) + n(- M )y} i+i = £ ( 2 J J(26-l)!![2( x f d - x L)] b (-2^r 



2/< 



Int(j/2) / \ fe 

cfy (2iy)^e ( ^ + ^ )2/(x ^- x ^ ) (31) 



These are the moments of a Gaussian distribution centered at ivj. Notice that \ud ~ Xud > 0, which follows from the 
explicit expression for the one- loop result (see Eq. §TU\i ). Since 

Imn(/ 1 ) = i(d( l u) + n(-/<)), (32) 
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the distribution of the imaginary part of the quark number is given by 

pi+i(y) = 1 e (.iy+»i) 2 /(xl d -x° d ) (33) 

ynxud - Xud) 

Note that the distribution takes complex values, as could be expected because of the phase of the fermion determinant. 
Moreover, we have that P^^ (y) is properly normalized, that the expectation value of the imaginary part of the quark 
number equals i times the isospin number in the phase quenched theory 



(Im[n]) 1+1 = / dyyP^ n] {y)=iv I , (34) 

J — oo 

and finally, that the average of the square of the imaginary part of the quark number is given by 

((Im[n]) 2 ) 1+1 = / dy y 2 P^ n] {y) = -uj + -( X L (35) 



In contrast to the real part, the width of the fluctuations of the imaginary part of n, (x^d ~ Xud) 1 ^ 2 > remains non zero 
for [i — > 0. 



C. The distribution of the quark number 

We finally turn to the distribution of the full quark number defined by 

^ +1 (^y) = ^^-^(«(M)-«(-^)))^y + 4( n ^) + «(-M)))) i ^ (36) 

within 1-loop chiral perturbation theory for fi < m^jl. In this section we show that this distribution factorizes into 
the distribution of the real and imaginary part of the quark number. 

Factorization occurs if the moments factorize. This can be easily seen by writing the probability distribution as an 
integral over the characteristic function, i.e. 

Pl +1 { X ,y) = J dk x dk y e -i^ e -iyk y < / e i^(n(^)-n(- M )) e ^W/ 1 )+«(- M ))\ ) ^^ (37) 

If the moments of the real and imaginary parts of the quark number factorize, the expectation values of the exponents 
in between the brackets will factorize. 

The relevant moments follow from a replicated generating function as follows 

((n(ji) + n(-f,)Y(n(p) - n(-,i)) fe ) 1+1 (38) 
= l im 1 d^ ■ ■ ■ d, j+k ( TT det(D + W 7o + m) ni det(D - W7b + m)"« TT detp + W7 ° + m) "' 



ij+k \f = \ jjv^ det(D - fxao + my 



1+1 



When we compute this replicated generating function in one-loop chiral perturbation theory a big simplification takes 
place: The contributions from Goldstone particles with one quark from the first two determinants and the other from 
one of the two determinants in the ratio exactly cancel. The reason is that the one-loop contribution of such mixed 
fermionic Goldstone particles occurs with the same combinatorial factor as the mixed bosonic Goldstone particles but 
with the opposite sign. This was first observed in Q where it was formulated as the absence of correlations between 
the phase factor and the magnitude of the fermion determinant to one loop order in chiral perturbation theory. Hence 
the moments of the real and imaginary part of the quark number factorize. 
The probability distribution therefore factorizes as 

Pi +1 (x,y) = P^(x)P^(y) 

ze -t*-vi? /(xL+x* d ) e (*y+vi) 2 /(xia-xS,) . (39) 

B \2 



^(xl d Y ix* d ) 



9 



To check this first main result let us first note that since the distributions of the real and imaginary parts of the quark 
number are normalized, also the product is normalized. Moreover, the expectation value of the quark number is zero 



Wi+i 



dxdy (x + iy)P„ 1+1 (x,y) = dx xPgfa (x) + i dy yP^ n] (y) = v T + iiui = 0. (40) 



We see that the total quark number (which necessarily is zero in chiral perturbation theory) is obtained only after a 
detailed cancellation between the contribution from the real part and the imaginary part. Such a detailed cancellation 
also occurs when we compute the average of n 2 



(n 2 )i+i = J dxdy (x + iy) 2 pi +1 (x,y) 

= Jdx x 2 P^{x) ~ jdy y 2 P^ n] (y) + 2 t Jdx xPg^x) J dy yP^ n] (y) 

= "i + \ixld + Xud) - (-"i + \(xid - Xud)) + Xviivi = Xud- (41) 

Note that, even though the distribution of the baryon number depends on the isospin density and off-diagonal 
susceptibility, these quantities drop out when evaluating the moments of the quark number operator. 

Finally, we note that since x u d + Xud ~~ * f° r A* ~~ * the quark number distribution becomes localized on the 
imaginary axis for /x = 0. This is in perfect agreement with the fact that the quark number operator is anti-hcrmitian 
for fi = cf. Eq. . 



IV. THE PQ DISTRIBUTION OF THE QUARK NUMBER OPERATOR 

In this section we give the result for the partially quenched distribution of the quark number. To derive this, all we 
need to notice is how the terms that mix the replica quarks and physical quarks contribute to the generating function. 

-Reai part: For the real part we consider Eq. (pt?|) where the expectation value is now taken in the 1 + 1* theory. 
Since Go (/i; , /i) is invariant under a change of the sign of both chemical potentials we get exactly the same mixing 
between the replica sector and the physical sector as before. Hence the final answer for the distribution of the real 
part of n is again the same 

^] ^ 
Imaginary part: The generating function for the imaginary part (|30p changes when we consider the partially 
quenched case. This time the two physical flavors make up an absolute square of the fermion determinant while the 
replica flavors makes up the phase factor. As we have seen before there are no correlations between these two factors 
within one-loop chiral perturbation theory. The generating function for the imaginary part of the quark number in 
the partially quenched case is thus given by Eq. (|30[) but without the single sum which mixes the replica and the 
physical sector. Hence there are no linear terms after differentiation and the final result is obtained from Eq. (|33p by 
setting Vj = 

R p Q(y) = 1 P -y 2 /ixL-x° d ) (43) 

\l<Xud-Xud) 

As it should (since we take vev's in the 1 + 1* theory) both distributions are real and positive. The full distribution 
is again the product of these two since the factorization only involves replicated flavors. Thus we find 

P^(x,y) = P^ n] (x)P^ [n] (y) 
1 



<x- V1 ) 2 /{ X I ui +Xud) p -v' 2 /{xU-Xud) 



(44) 



As a cross check we sec that 



^{xLY - (xSd) 2 

(n p Q) 1+1 ,=^ (45) 



and 



((n PQ ) 2 h + i*=X B ad + vl (46) 

are in agreement with table Q] 

In section [VT] we make use of these results when discussing the problems faced by numerical lattice QCD at fi ^ 0. 
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V. THE DISTRIBUTION OF THE CHIRAL CONDENSATE 



In this section we derive the distribution of a and study how the chiral condensate (a) builds up. As was the case 
for the quark number, the operator a is not hermitian 

a{uY = ( Tr ^ = Tr = a(-fi). (47) 

\ D + /i7o + m J D - fx'jo + m 

Therefore we derive the distribution in the complex a plane. We start by computing the distribution of the real and 
of the imaginary part separately. The distribution of the full chiral condensate then follows as the product of the two. 
This is precisely the same which happened for the baryon density. In fact, the derivation for a is almost identical to 
the one for the quark number. 

To see that the derivation is analogous to that for the baryon number let us start with the distribution of Re[er] 
and derive the generating function. We evaluate the distribution 

P£J, (*) = (s(x- \{a{n) + i > C 48 ) 

for /i < m n /2 to one-loop order in chiral perturbation theory. 

The 5-function can be represented as an integral 



pi+i 

Re[<r] 



1 f / \ 

(a) = _ / dk e~ lxk ( e +** ) (49) 

J 2?r J_ oa \ 1 1+1 

and the exponential is expanded 



e 4 *W*>+-H0>) = E^((4(^)+^(-M))) ) , (50) 

1 + 1 3=0 ' I 1+1 

so that the distribution follows from the moments. 

The moments can again be expressed in terms of a replicated generating function. This time the replica index 
labels the masses m,\ , . . . , rrij , 

(a(ri+a(-[x)y) (51) 
/ i+i 

= lim d mi ■■■d mi ( TT det(D + fijo + mi)"' {D - /no + m») ni ) 

«i^o m ■ ■ ■ rij \ / 

J \i=i i i_| 

We now need to keep track of the mass dependence of the replicated generating function 

/ JJ det(U + ^70 + m,)"' (D - M7o + ) (52) 

\i=l / 1 + 1 

( 3 

= cxp 2J 2n i n m (G ! o(jU,M ! mi,m m ) + G (-M,/x,m;,m m )) 

\l<m=l 

+ ^2 2n i \ m i^ + ( g o(m, m i m i) + Go(-fJ>; A*, m, mi))} . 
i=i ) 

Compared to the generating function for Re[n] in Eq. (|120"|) there are two differences: 1) The replica index now labels 
the quark masses instead of the chemical potentials. 2) The term linear in the replica number contains the mean 
field value, E, of the chiral condensate. However, the entire structure of the generating function remains and the 
combinatorics associated with the differentiation is identical to that for the quark number. Only now the physical 
quantities appearing in the expressions arc the chiral condensates and the chiral susceptibilities. We get 

P^ix) = - 1 e-fr-s-'s'-^V^+xf), (53) 
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where we used the definitions in (pi 

Likewise the distribution for the imaginary part of a 

P^M (y) = (s(y + 4(cr( M ) - "(-P))))^ - ( 54 ) 

follows form that for the imaginary part of the quark number simply by re-identifying the physical quantities which 
appear in the final expression 

pWjQ,) = 1 e fa-^W)V(x|-xf), (55) 

Notice the absence of the mean field value of the chiral condensate. In the generating function the mean field 
contribution from the fcrmionic replicas cancels against the one from the replicated bosonic quarks. In agreement 
with the hermiticity property (|47|) we see that the width of vanishes for /j, — > while the width of the distribution 

of the real part remains nonzero. From the explicit expression for Go it can be easily shown that Xs — x$ > 0- 
Since the structure of the generating function is unchanged the full distribution, 



again factorizes 



Thus we have 



Pl +1 {x, y) = {5[x- -(a(p) + ) 6 + %-{a{ji) - <j(-^))jj^ , (56) 

Pl +1 {x,y) = P^ ] {x)P^ a] {y). (57) 

P 1+1 (x y) = - e -{x-s-5T, B -5T, 1 f/{x l s+xs ) e (*y-tz B +SS 1 ) 2 /( x ' s - x f) 

v^ix'sY (xf ) 2 

As a check of this second main result we compute the chiral condensate 

/oo 
dxdy (x + iy)P^ +1 (x,y) =£ + 2<5£ s , (59) 
-OO 

and the square of the condensate 

poo 

(* 2 )i+i= dxdy (x + iy) 2 pi +1 (x,y)= X § + (S + 2<5I] S ) 2 . (60) 



We see that the dependence on S 7 and Xs nas canceled and the results obtained in table U are reproduced. The 
cancellations take place in exactly the same manner as for the baryon density. We conclude that in both cases the 
fi dependence of the fluctuations are induced through a coupling to the isospin charge of the pions. These strong 
fluctuations in the real part and in the imaginary part combine and leave the physical observable independent of the 
chemical potential. 

Finally, we give also the partially quenched distribution of the chiral condensate 

pPQ(x,y) = 1 _ p -(x-s-ra*-ra')'/(xS+xf),-v a /(xS-xf) (61) 

^VW-(xf) 2 

The partially quenched expectation values from table U 

(<j) PQ =S + <5E B + 5E / , (62) 

and 

<^ 2 >PQ=xf + (Z + SX B + SX 1 ) 2 (63) 

follow from Gaussian integrations over the distribution function. 

Wc will use both the partially quenched results and the full results to explain how to deal with some of the numerical 
problems encountered in lattice QCD at nonzero chemical potential. 
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VI. NUMERICAL LATTICE QCD AT NON ZERO fi 

Above we have derived the full and partially quenched distributions of the quark number operator for /i < m T /2. 
Here we discuss how these results can be of use when measuring the quark number operator by numerical lattice 
QCD. 

The main problem encountered in numerical lattice QCD at nonzero chemical potential is that operators may 
acquire their expectation values by virtue of extremely delicate cancellations caused by the complex valued fcrmion 
determinant in the path integral. To illustrate the problem let us look at the distribution of the quark number 
operator, Eq. ([39)) . Notice that the distribution itself takes complex values, i.e. it is not a probability measure. This 
feature is also shared by the unquenched eigenvalue density of the Dirac operator [32l . |33| and the distribution of the 
phase of the fermion determinant These complex oscillations are all important in order to separate the physics 
of nonzero baryon chemical potential [H| from that at isospin chemical potential. Common to the three examples 
mentioned above is that the amplitude of the complex oscillations is exponentially large in the volume. Moreover, in 
all three cases, one must integrate over on the order of V periods of the oscillations in order to obtain a reliable value 
of the baryon density or chiral condensate. Let us illustrate this explicitly using the results derived above. 

The unquenched distribution of the quark number, P^ +1 {x,y), takes complex values, and, as we have seen in 
Eq. (|40|) . both the real part and the imaginary part contribute to the baryon density. The magnitudes are equal, fj, 
but the signs are opposite so that the total quark number vanishes. Since the complex oscillations are associated with 
the imaginary part, y, of the quark number, n = x + iy, let us ask: How large should y ma x be in order that 

ji+i 



dyiyP^ n] {y)~- Vl . (64) 
The answer is: 

vliax - v\ > xld ~ Xud (65) 

Since vi ~ V and x u d ~ Xud ~ ^ we find that y m ax has to be only slightly larger than vj. As the period of the 
oscillations is of order unity we conclude that we have to include on the order of V oscillations in the integral. This can 
be a hard task to control numerically unless we know the analytical form of the distributions. The results presented 
here give this form to leading order. There will be corrections to this form in order to induce a nonzero baryon 
number, but one should still expect large cancellations between contributions from the real part and the imaginary 
part. 

Not only can results from chiral perturbation theory help to understand the detailed cancellations occurring in 
the integral, they can also give hints on which part of the integrand one needs to sample numerically in lattice 
QCD: Imagine that we have generated an ensemble of configurations for the phase quenched weight. Then for the 
configurations in our simulation ensemble, the quark number n is distributed according to P^®(x,y) as given in 
Eq. ([44]) . Within the width of this distribution only on the order of vV oscillations takes place. The reweighting from 
the simulation ensemble to the full theory (target ensemble) therefore has to lift the far tail of the phase quenched 
distribution. (For illustration, see figure [1]). 

Even though, the sign problem makes it hard to approach the thermodynamic limit V — > oo, it is worth to keep 
in mind that in a given numerical simulation we work with a finite volume where the sign problem may be tractable 
depending on the value of T and ji. We will address this issue in [35[ ■ 

In the Complex Langevin [22, HH, H3, [H, [2(| approach to the QCD sign problem the real and the imaginary part of 
the quark number are themselves complex. In particular, this means that the imaginary part, y, is not constrained to 
the real axis. Rather as we now show, y will fluctuate parallel to the real y axis but shifted into the complex y plane 
by iv\. To show this we start from the complex action, 

S = - log[P^](i/)] = -(iy + vi? I {xld ~ X B ud), (66) 

where we made use of the result in Eq. (|33|) . The flow equations for y = a + ib are given by (the step size is denoted 
by e) 



a-n+i — a n — eRc 



dS_ 

dy 



y=a n -\-ibn 



a n - e j 2 °" B + V«7„ (67) 

X u d X u d 
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FIG. 1: Left: The distribution of the imaginary part of the quark number in the partially quenched simulation ensemble. 
Right: The contribution to the quark number in the 1+1 target theory from the imaginary part of n as a function of y m ax, 
see Eq. (|64|l . Only when y max is larger than Vj do we find that the contribution from the imaginary part cancels the one from 
the real part and leaves (n)i+i =0. If vi ~ V is much larger than the width of the PQ distribution \J \ ud — Xud ~ V* j then 
there is a serious overlap problem. For the plot we have chosen vi = 20 and Xud ~ Xud = 10; an d we already see that, unless 
one also samples the far end of the tail, the configurations generated in the PQ simulation ensemble do not have the values of 
Im[n] required for the cancellation of the baryon number. 



and 

h h T \ dS ] 

b n+ i = b n - elm — 

. a V 1 y=a n +ib„ 

2(b n - vi) 

= o n - e j - g . (68) 

X u d X u d 

Note that a and b decouple: a fluctuates about zero while b quickly moves to i>i and stays there, since there is no 
noise 77 to kick it around. The complex Langevin algorithm, therefore, will essentially replace 

00 

<[!](?/) (69) 

-OO 

by 

da {a + wi)e" al/{x '^- x ^\ (70) 



where a max is the width of the region sampled by the algorithm. Since the flow equation for a is just that of the 
ordinary real Gaussian, cxp(— a 2 /{x u d ~ Xud))' anc ^ completely decoupled from the imaginary part, the complex 
Langevin algorithm should not have any problems in sampling this. In other words, a max is a good deal larger than 



Xud ~ Xud an( l t nc integral is close to ivi as desired. Again there will be corrections to the one-loop results derived 
here and complex Langevin must be able to take these into account correctly in order to obtain the correct average 
quark number. However, as long as the average quark number is much smaller than vi the fluctuations in the complex 
y plane should be expected to take place in the neighborhood of ivi. A similar example of a possible usefulness of 
complex Langevin was presented in [36j. 

For purely imaginary values of the chemical potential the sign problem is absent, and the quark number operator 
is imaginary. In |Appcndix A| we work out the distribution of the quark number for imaginary values of the chemical 
potential. The result is a Gaussian centered at zero. It is certainly true that simulations at imaginary /x (37l . I38L l39l | 
are easier - the reason is the positivity of the measure. This suggests that simulations at imaginary fj. are as easy 
as simulations at [i = 0. Analytical and numerical studies show that the analytic continuation is well under control 
at least for fi/T < 1 [3?], HH, |23|, and it would be interesting to interpret this result in the light of the distributions 
discussed here. Analytic continuation to real fi has also been studied in models without the sign problem (40| . 
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VII. DISTRIBUTIONS FOR fj, > 



The chiral expansion of the generating functions relevant for the baryon distribution take a different form when 
/i > mjr/2. In this domain the chemical potential induces a Bose condensate with a nonzero isospin number even at 
the mean field level. The derivation presented above is therefore not valid for p > m w /2. In this section we discuss 
the distributions of the quark number and the chiral condensate for p > m^/2. 



A. The fluctuations of the chiral condensate (p > m^/2') 

For p > 771^/2 the quark mass is inside the spectrum of D + pjo [HI, [U [42|, H^|. Since an eigenvalue of D + pjo 
can come very close to the quark mass the fluctuations of the chiral condensate are much larger when p > m w /2. In 
order to quantify the fluctuations of the chiral condensate let us consider the moments of the real and the imaginary 
part of a. To start we consider the quenched case. 

The odd moments of Im[er]: Because the quenched weight does not depend on the sign of p the odd moments vanish 

({a - a*) 2p+1 ) =0. (71) 

It follows that the quenched distribution of the imaginary part of the chiral condensate must be an even function; 
that is: -P I ^j^.~°' ) (y) is symmetric in y. 

The even moments oflm[a]: The even moments are nonzero. In fact, as we now show, they are divergent. To see 
this let us first consider the second moment of the imaginary part 

((a-a*) 2 )=2(a 2 )-2(aa*). (72) 

If we express a in terms of the eigenvalues, Zk, of D + /170 

* = r- ( 73 ) 

k 

we see that (crer*) includes an absolute squared pole, which gives rise to a logarithmic singularity. So we have 

{(a-a*) 2 ) = -2(a<T*) + O{e ). (74) 

Obviously, the logarithmic singularity is only present when the quark mass is inside the support eigenvalue density. 
In fact, it can be shown explicitly within chiral Random Matrix Theory that the divergent part of the second moment 
is proportional to the eigenvalue density evaluated at the quark mass 



<(^*) 2 > = - 2 (E i z i 2 ) = - 2 / fzpiz^T-tr-* 

k l Zfc + m l / Jc c \z + my 



VT, 2 

= 4nlog(e)p Nf = (z = m,z* = m) = log(e)6>(>| -m n /2)-^. (75) 

Here, C e is the complex plane excised by a sphere of radius e centered at — m. The quenched eigenvalue density is the 
leading order result from chiral perturbation theory [42j. We conclude that the quenched second moment of Im[ir] is 
logarithmically divergent if the quark mass is inside the support of the eigenvalue density, i.e. if [i > m n /2. 

The higher even moments of a have more severe divergences. The most divergent term is the one with the same 
powers of a and a*, so that, 



((a-O^^l-l^K)'). (76) 
plpl 



In order to understand this better let us work out the details explicitly for p = 2. Again we keep only the leading 
divergence 

((a -a*) 4 ) = 6((aa*) 2 ) 

= «(j:— — - 1 / V (77) 

\ z — ' z a + m Zh + m z* + m zf + m / 

\g,h,j,k y J I 
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The most singular terms are the diagonal terms with g = h = j = k which lead to a 1/e 2 singularity which dominates 
the integral over the spectral density. 

Because the most divergent term is also given by the diagonal part of the sum for p > 2, it is always proportional 
to the eigenvalue density. We thus obtain the general relation (for p = 2, 3, 4, . . .) 

We conclude that P^f^ (y) must be symmetric in y in order that the odd moments vanish, and it must have a 

l/|y| 3 tail to reproduce the observed singularities. When the quark mass is outside the spectrum of the Dirac operator 
the same argument applies, but now the spectral density at z — m is exponentially suppressed, and the singular terms 
vanish in the thermodynamic limit. While the singularities tell us about the tail of the distribution, the leading 
divergent behavior of the moments docs not contain sufficient information to obtain Pyj" (j) for smaller y. We 

have verified the l/|y| 3 tail of (y) by a numerical simulation of quenched chiral Random Matrix Theory. 

So far we have considered the quenched case. As a first step toward the fully unquenched theory let us consider the 
phase quenched theory. In this case the absolute square of the fermion determinant in the measure, Y[j \ z j + m2 \ 2 ■> 
shifts the singularities. We now have 

((a-a*) 2 ) 1+1 , ~e°, (79) 



<(a-a*) 4 ) 1+1 , ~log(e) (80) 

and 

(81) 

for p > 2. This implies that the far tail of the distribution of Im[a] drops of like 1 / 1 2/ 1 5 - The distribution is again even 
since the odd moments still vanish because the absolute square of the fermion determinant does not depend on the 
sign of the chemical potential. 

Finally, let us briefly look at the unquenched case. Let us for simplicity take Nf = 1. With a single determinant in 
the measure the odd moments no longer vanishes. Moreover, the odd moments can now also diverge. For example, 

((a - <r*) 3 )N f =i = -(a 2 a*) Nf=1 ~ log(e). (82) 

Again the log(e) singularity is only present when the quark mass is inside the support of the spectral density, i.e. for 
/i > m n /2. In general we have 

((cr-O^W-^iy (83) 

for p > 1. The diverging odd moments show that the unquenched distribution of Im[cr] is not an even function of y. 
We expect that the unquenched distribution of Im[er] takes complex values. 

Unfortunately, due to the divergences, we have not been able to derive the full distribution of the imaginary part 
of the chiral condensate for /i > m n /2. 

For the real part of a the main difference is that the odd moments are nonzero. The leading divergent part of the 
even moments is the same (up to a sign) as for the imaginary part. 



B. Fluctuations of the quark number (/i > /2) 



The general arguments which gave us information about the distribution of the chiral condensate also apply to the 
quark number for p > m^/2. For p > m^/2 we have that \x is inside the support of the spectrum of 7o(-D + m) [iil ]. 
and this again leads to enhanced fluctuations. A rerun of the general arguments presented above show that also the 
quenched distribution of the real and of the imaginary part of the baryon number has an inverse cubic tail. Rather 
than repeating the arguments let us instead discuss the cancellations which insure that the total baryon number 
remains zero in the unquenched case (as must be true in chiral perturbation theory) . 
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Let us first consider the real and imaginary part of the quark number density for the quenched or phase quenched 
case. The imaginary part of the quark number is given by 

2i Im[n] = Tr ]° + Tr — g . (84) 

Because the spectrum of 70 (-D + m) is reflection symmetric about the imaginary axis, and the average spectrum for 
the quenched and phase quenched theory is reflection symmetric both about the real and the imaginary axis, we find 
that also for fi > 

(Im[n]) Nf=0 = (Im[n]) 1+1 , = 0. (85) 
The real part of the quark number is given by 

2Re[n] = Tr— 2° Tr— . (86) 

7o(-D + m)+/i 7o(Z? + to)-^ 

This is the isospin density with a nonvanishing expectation value in the quenched and phase quenched theory for 
fi > because of pion condensation. At mean field level in Chiral Perturbation Theory the isospin density is [47| 

^ = (Rc[n]) 1+1 , = 2/ iF 2 [l-(^) 2 ]. (87) 

We now consider the unquenched case with Nf = 2. In the supersymmctric formalism, the generating function for 
Re[n] is given by 

/ det( 7o (£> + to) + ^ 2 )det(7 (D + m) - /i 2 )det 2 (7 (£> + m) + m) \ 
A{(11 ' fX2) \ det(7o(£> + m) + //i)det(7o(i) + m)-Mi) / 1 ' 



with 

(89) 



2(Re[n])i+i = -j— log Z(jjl 1} /x 2 ) 
dfi 2 



The generating function undergoes a phase transition to a pion condensed phase at fj, = m w /2. It can be interpreted 
as the average phase factor at /ii for the phase quenched theory at fj.2- Such averages where studied in [45[. At the 
mean field level the generating function factorizes as 

7{ x = (det(7o(£> + m) + fi 2 )dct(-yo(D + m) - i u 2 )dct 2 (7o(£' + m) + m)) 
(flufl2) - (det(7o(D + m) + /xi)det(7o(£> + m)-Mi)) ' 

The imaginary part of the quark number for two flavors can also be obtained from the generating function 

d 



2i(lm[n]) 1+1 

a Hi 



logZ( Ml , M2 ) (91) 



with 



/det(7o(-D + m) +/i 1 )det(7o(D + m) +/x 2 )det(7o(-D + m) -/i 2 )\ , nn , 

a\HuH2) = \ ttt — 777-; \ s /• 

\ det(7o(-D + m) - /ii) / 

The partition function (|9"2")l is a phase- quenched average phase factor which was studied in (38[. At the mean field 
level this partition function factorizes as 

= (dct 2 (7o(£) + to) + m) det(7o(£> + m) + fi 2 ) dct(j (D + m) - M2)) 
lMl,M2j " (det(7o(D + m)-Mi)det('yd(D + m)+/ii)> ' 

We observe that the real and imaginary part of the quark number are obtained from the same generating function. 
The mean field analysis of partition function (f9T)]) was outlined in (4|| . It is determined by the action 

5= -F 2 Tt[U,B][U- 1 ,B]-TrME{U + U- 1 ) (94) 
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with the baryon matrix and the mass matrix given by 

B = diag(/ii, jUi, /i2, —^2), M = diag(m, m, m, m). (95) 
The ansatz for the mean field can be written as [46| 

U = R^(l3 k i)R ik {a k )R k i(f3 kl ) (96) 

with k = 1, 2, 3 and I ^= k. Here, R pq (a) is a rotation in the pq plane by angle a. There is one important difference 
with [46| though. Because the quark masses are equal and the chemical potential are put equal after differentiation, 
the integration over (3 k cannot be done by a saddle point approximation but has to be performed exactly. 

In total there are 6 different saddle points. For /ii = \ii the action of each of the saddle point is the same and the 
dependence on f} k cancels. Therefore, the integral over (3 k has to be performed exactly for [i\ ^ ^2, whereas the a k 
are determined by the saddle point equation. 

For the real part of the quark number we obtain: 

2(Re[n])i +1 = -(12 + 4)c^F 2 [cos 2 a - 1], (97) 

where the first term (i.e. 12) originates from the rotation matrices that mix — ^2,^21/^1 arL d the second term (i.e. 4) 
from the rotation matrices that mix — ^jMliMl- The real constant c is a normalization factor. The solution of the 
saddle point equation is given by (4?j] 

77Z 2 

cosa=-^. (98) 
4/x 2 

With this we have 

2 

Tfi 

(Re[n]) 1+1 =8 CM F 2 [l-( i ^) 2 ], (99) 

which is proportional to the isospin density in the Bose condensed phase of the phase quenched theory at mean field 
level,cf. (|87j). 

For the imaginary part of the quark number we obtain contributions both from the numerator and the denominator 
of the generating function. The contribution that originates from the numerator is given by 

(4 + 4)/iF 2 [cos 2 a- 1]. (100) 

where again the first and second term correspond to rotation matrices that mix —^2, ^2, Mi an d that mix — /12, fJ-i, ^1 
in this order. The contribution from the denominator is given by 

- 24c,uF 2 [cos 2 a- 1]. (101) 

Since the real and imaginary parts are obtained from the same generating function, the constant c in the expressions 
is the same. For the imaginary part of the quark number we thus obtain 

2i(Im[n])i+i = 16c^F 2 [cos 2 a - 1]. (102) 

As must be true in chiral perturbation theory the sum of the real and imaginary part of the quark number vanishes. 
Here we have shown this by explicitly computing both contributions. As for /i < m n /2 we have found that the two 
terms arc proportional to the isospin density in the phase quenched theory. For fi > m n /2, however, the isospin 
density is far greater due to Bose condensation of pions. 



VIII. CONCLUSIONS 



The distribution of the quark number operator and the chiral condensate for Euclidean QCD at nonzero chemical 
potential has been derived to leading order in chiral perturbation theory As the two operators take on complex values, 
the distributions are over the complex plane. Moreover, because of the phase factor of the fcrmion determinant, the 
distributions are not real and positive. We have shown how the complex oscillations of the unquenched distributions 
lead to large cancellations when evaluating the baryon density and chiral condensate. These cancellations give a direct 
insight into the problems faced by numerical lattice QCD at nonzero chemical potential. Of course the net contribution 
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from pions to the average baryon density is zero within chiral perturbation theory. Nevertheless, pions contribute 
manifestly to the distribution of the quark number and the chiral condensate and hence to the noise produced in 
numerical lattice gauge simulations at nonzero chemical potential. 

Most of the results were derived for /i < to one- loop order in chiral perturbation theory. Then the distributions 

of the chiral condensate and the quark number take a Gaussian form as one might expect from the central limit 
theorem. The unquenched distribution of the imaginary part of these observables takes on complex values. For 
example for the imaginary part of the quark number, this distribution is a Gaussian that is shifted in the imaginary 
direction by an amount proportional to the volume (the isospin number in the phase quenched theory when evaluated 
at the same value of the chemical potential) . This is of course not possible in an ordinary reweighting scheme and the 
final results for the baryon density instead relies on detailed cancellations after taking into account the far tail of the 
distribution. However, within the Complex Langevin method the real and imaginary parts of the baryon density are 
complexified and the imaginary part could fluctuate around i times a quantity of order the volume. 

For [i > TOtt/2 the distributions of the real and the imaginary part of the quark number and chiral condensate 
develop a power law tail. This extreme enhancement of the fluctuations is a direct consequence of the quark mass 
being inside the spectral support of the Dirac operator, D + /i7o, and the chemical potential being inside the support 
of 7o(-D + to). 

The analytical results provided here may help in interpreting state of the art lattice simulations at non zero chemical 
potential. They can be used as benchmark for attempts to link the fluctuations of the baryon number to the presence 
of the tricritical point. Since we have understood fluctuations produced the pions this may help in optimizing and 
developing numerical approaches to reduce this source of noise. 
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Appendix A. IMAGINARY CHEMICAL POTENTIAL. 



In this appendix we give the result for the distribution of the quark number operator when evaluated at purely 
imaginary values of the chemical potential. The fermion determinant is real for imaginary chemical potential, and 
numerical simulations are possible for an even number of flavors (37l . l3Sl l39j . In this case the quark number operator 
is anti-hermitian even when i[i is nonzero 

n(«»* = ( Tr 25 ^ = -Tr ^ = - n {i^). (103) 

\ D + i/i7o +mj D + ifjrfo + m 

The distribution is therefore one-dimensional 



^(y) = (S(y + in(i^)) 1+ltill . (104) 



The derivation of the distribution for fi < to one-loop order in chiral perturbation theory is simpler than for 

real \x since now there are no terms n{—i\i) which can couple to the isospin charge of the pions. We find 

CM = 7= eV/li ' ( 105 ) 

As must be true within chiral perturbation theory we find that 

(»M>i+i, v = 0- (106) 

Note also that the width of the distribution is independent of [i. As Xud 1S extensive the distribution of the quark 
density becomes a (5-function at the origin in the thermodynamic limit. 
We also get 

(nM 2 ) l+lii)1 - \Xud\- (107) 

One can of course also consider the distributions of n(i^i) ± n(—ifi) which will be the analytic continuations of the 
distributions of the analytic part of the real and imaginary parts of n. 
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